Lp-THEORY FOR A CAHN-HILLIARD-GURTIN SYSTEM 

"^" I MATHIAS WILKE 

o 

J^ ■ Abstract. In this paper we study a generalized Cahn-Hilliard equation which 

was proposed by Gurtin [8]. We prove the existence and uniqueness of a local- 
in-time solution for a quasilinear version, that is, if the coefficients depend on 
j<^». , the solution and its gradient. Moreover we show that local solutions to the 

^v/l ' corresponding seinilinear problem exist globally as long as the physical potential 

satisfies certain growth conditions. Finally we study the long-time behaviour 
of the solutions and show that each solution converges to a equilibrium as time 
P^ , tends to infinity. 

• ' 1. Introduction 

"t^ . We start with the derivation of the classical Cahn-Hilliard equation. Consider the 

free energy functional of the form 



- 1 — \ 



/, O"^ 



(1.1) j-(^)^y^ (_|vv.|^-f$(v,)l dx, 

where fi is a bounded, open and connected subset of M" with boundary F := dil G 
C^. We assume that the order parameter V' is a conserved quantity. The according 
^^ , conservation law reads 

^^ . (1.2) dt^ + divj=0, 

CO . where j is a vector field representing the phase flux of the order parameter. The 

^— ^ ' next step is to combine the two quantities j and /i. Similar to Fourier's law in the 

^ ^ . derivation of the heat equation one typically assumes that j is given by 

(1.3) J = -Vm, 

a postulated relation. Finally we have to derive an equation for fi. The chemical 
potential /i is given by the variational derivative of J^, i.e. 

<^ • SF 

^ = — = -A^ + $'(V'). 

If F is of the form (jl.ll) this yields the classical Cahn-Hilliard equation. 

In the early nineties GuRTiN [5j proposed a generalized Cahn-Hilliard equation, which 

is based on the following objections: 

• Fundamental physical laws should account for the work associated with each 
operative kinematical process; 

• There is no clear separation of the balance law (|1.2p and the constitutive 
equation (|1.3p : 
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• Forces that are associated with microscopic configurations of atoms are not 
considered in the derivation of the classical Cahn-Hilliard equation. 
According to Gurtin there should exist so called 'microforces' whose work accompanies 
changes in the order parameter tp. The microforce system is characterized by the 
microstress ^ e M" and scalar quantities tt and 7 which represent internal and external 
microforces, respectively. The main assumption in ^8 is that ^, n and 7 satisfy the 
(local) microforce balance 

(1.4) div^ + 7r + 7 = 0, 

which can be motivated from a static point of view, see [5] for more details. In a next 
step we want to derive constitutive equations, which relate the quantities j, the flux of 
the order parameter, ^ and tt to the fields tp and fi. The technique used in f8| for this 
derivation is based on the balance equation ()1.4|) and a (local) dissipation inequality, 
which is a direct consequence of the first and the second law of thermodynamics, that 
is, the energy balance 

e dx ^ - q-vda+ I r dx + W {Vl) + M{Vl) , 



di Jn Jdn Jn 

and 

— S dx>- - ■ v da + - dx, 
"s Jn Jdn f Jn^ 

cf . [5J Appendix A] . The second law of thermodynamics is also known as the Clausius- 

Duhem inequality. Here e is the internal energy, S is the entropy, is the absolute 

temperature, q is the heat flux, r is the heat supply, W(n) is the rate of working 

on ft of all forces exterior to i} and A4{^) is the rate at which energy is added to 

fi by mass transport. Let F be the free energy density, depending on the vector 

z = {tjjy'Vijj., fi,\7 fijdt^p)- Then the second law of thermodynamics (in its mechanical 

version as considered by Gurtin |Sj) reads 

— / F{z) dx < — l^j{z) ■ v dcr + I ^ • vdiip da + fim dx + jdtip dx, 
dt Jn Jdn Jdn Jn Jn 

with m being the external mass supply. Making use of Green's formula, we obtain 



— / F{z) dx<~ {\'^I■J{z) + ^idivj) dx 
"t Jn Jn 



+ / {div^dti^ + ^ ■ Vdt^j) dx + fim dx + jdt'tjj dx. 
Jn Jn Jn 



in presence of external mass supply m, ()1.2p will be modified to 

(1.5) dtip + div j = m. 

In view of (|1.4I) and (|1.5p we obtain the dissipation inequality 

— / F{z) dx < {pdtip — j ■ VfJ- — Trdtijj + £, ■ ydttp) dx. 
dt Jn Jn 

This in turn yields the following local dissipation inequality 

dtF{z) < ^idti! - J • Vm - Tidti^ + C • VStV, 
for all fields ip and /x, this means, we have 

(1.6) (a^F + vr - /i)^ + {d^^F - • Vt/^ + 9^F/i + d^f^FWfi + d^F^, + V^i • j < 0, 
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where u = dtu and u — dfu for a smooth function u. This local inequality needs to 
be satisfied for all smooth fields ^ and /i. Hence we have necessarily 

F{z) ^ F{i,,V^) and ^IV-, V^) = 5v^F(V^, Vi/^) 

and there remains the inequality 

{d^F + TT- ^i)^j + V/i • j < 

whose general solution is given by (cf. (H Appendix B]) 

9,0_F + n ~ fi — —/?■(/' ^ c • V/i and j = —aip — B\7 ii, 

with constitutive moduli (3{z) (scalar), a{z),c{z) (vectors), B(z) (matrix) and the 
constraint that the matrix 

B 



(1.7) 



is positive semidefinite. For convenience we assume that /3 is constant and a, c and 
B do only depend on x instead of z, whence we deal with an approximation of the 
constitutive moduli l3{z),a{z),c{z),B{z). In particular, if the free energy density F 
is given by F{ip,Wtp) — i|V'(/'P + *&(''/') '^^ obtain the following semilinear Cahn- 
Hilliard- Gurtin equations. 

dt-ip - div(SV/i) - div(aatV') = /, t>Q, xeVl, 
^'' /i-c- V/1 + AV'-/39tV'-$'(^) =5, i>0, xel7, 

where fi C R" is open and bounded with boundary F := dO, E C^. We want to 
emphasize that for the special case B = I,a = c = and /3 = 0, we obtain the 
classical Cahn-Hilliard equation or the viscous Cahn-Hilliard equation if /3 > 0. 
Let us point out that we will also deal with a quasilinear version of (jl.81) in Section 
5. To be precise, we will consider the system 

dtip - div(6(2:, ■0, ^ipWlJ-) - div(a(2;, ip, V^)dt^) = /, i > 0, a; e fi, 
(1 9) 
^ ' ' PL- c(x, -0, VV') • V^i + A^/i - Pdtip - $'(-0) =5, i > 0, a; e f^. 



In this paper, we are interested in solutions of (jl.Sp and ()1.9p subject to the Neumann 
boundary conditions for [ip^p), having optimal Lp-regularity in the sense 

i, e Hl{j-Hl{n))nLp{j-Hl{n)), 

and 

MeLp(J;Hp2(r!)), 

for given functions / € Lp{J;Lp{n)) and g £ Lp{J; Hp{^)), where J = [0,T]. We 
will always use the following assumptions for the semilinear problem 

• a,ce CH^)", 

• diva(a;) = divc(a::) = in 51, 

• {a{x)\i^{x)) — {c{x)\i'{x)) = on dV,, 

• I3>0, B = bl, with b e C^m, 

• there is a constant e > 0, such that the estimate 

I3zl + {a + c\zi)zo + (Bzilzi) > e(z^ + jzip) 
is valid for all {zo,zi) G R x M" and all x £ ^. 
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In Section 2, where we consider il = M", we allow for general, positive definite matrices 

B. It is also possible to consider those matrices in all other sections but for the sake 

of convenience we restrict ourself to the case B = bl. Actually this allows to draw 

back the problem in the half space R" to the whole space R" by means of reflection 

methods. 

Results on existence and uniqueness can be found e.g. in the papers of Bonfoh & 

MiRANVILLE 3°, MiRANVILLE [lO], [ll], MiRANVILLE & PlETRUS 16], MiRANVILLE, 

PiETRUS & Rakotoson [T2] and MiRANVILLE & Zelik 14j. In any of these papers 
the authors use a variational approach and energy estimates to obtain global well- 
posedness in an L2-setting, with periodic boundary conditions for a cuboid in R"^. 
The qualitative behavior of solutions of the Cahn-Hilliard-Gurtin equation has been 
investigated in f^, [12] and [13]. In [3] and [12] the authors proved the existence 
of finite dimensional attractors, whereas MiRANVILLE & Rougirel 13} showed that 
each solution converges to a steady state, again with the help of the Lojasiewicz-Simon 
inequality. One assumption of Miranville & Rougirel [13] is that the norms \a\, 
\c\ and |i3 — /| are bounded by a possibly small constant. In the present paper we will 
give an alternative proof for the relative compactness of the orbit {ip(t)}t>o in i/^ (^) 
with the help of semigroup theory and a priori estimates (see Proposition [TH]). 
The present paper is structured as follows. In Section 2 we deal with a corresponding 
linearized system to (|1.9p in the full space R" with constant coefficients. Section 3 
is devoted to the analysis of the linearized system with constant coefficients in the 
half space E". Making use of the optimal regularity results of Sections 2 and 3 
we apply the method of localization and some perturbation results in Section 4 to 
derive optimal Lp-regularity for the linearized Cahn-Hilliard-Gurtin equations (i.e. 
p.Sp with $' = 0) in an arbitrary bounded domain fl C R" with boundary dil € C^ . 
In Section 5 we prove the existence and uniqueness of a local-in-time solution of (jl.9p . 
For this purpose it is crucial to have the optimal Lp-regularity result from Section 
4 at our disposal. To the knowledge of the author there are no results on the local 
well-posedness of (|1.9p but only for the case where a, c and h depend solely on the 
order parameter -0, cf. Miranville [TF. In Section 6 we investigate the global well- 
posedness of the semilinear system (|1.8I) . The basic tools are a priori estimates and 
the Gagliardo-Nirenberg inequality. Finally, in Section 7, we show that each solution 
il){t) of ()1.8|) converges to a steady state in H2{^) as i — ?> oo. To this end we will use 
relative compactness results and the Lojasiewicz-Simon inequality. 

2. The Linear Cahn-Hilliard-Gurtin Problem in M" 

In this section we will solve the full space problem 

dtu - dbj{adtu) = div{B\7^) + f^ t > 0, a; G M", 

(2.1) n-c-V^i^ /3dtu - Au + g, t > 0, x e R", 

u(0) = uo, t = 0, xe R", 

where ,3 > 0, a,c e R" and B e R"""'. Note that the matrix (fLT]) is positive 
semidefinite if and only if 

/3z'^ + {a + c\zi)zQ + {Bzi\zi) > 

holds for all (zo,zi) € R x M" and all x G il. Here (-j-) denotes the usual scalar 
product in C" and the vector fields a, c as well as the matrix valued function B are 
assumed to be smooth. In the sequel we will use a slightly stronger assumption. 
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(H) There is a constant e > 0, such that 

Pzl + (a + c|zi)zo + {Bzi\zi) > e{zQ + \zif) 
is vahd for all {za,zi) e M x M" and all a; G fi. 

The following result is useful for the analysis of (j2.ip (see also [131 Lemma 5.1]). 
Proposition 2.1. Let (H) hold. Then 

mm -\{{a®c + c® a)m > em^ 

for all^eW. 

Proof. Hypothesis (H) reads 

Pzl + (d|zi)zo + {Bzi\zi) > e{zi + \zi\^), 
where d := a + c. Observe that the left side of this inequality can be rewritten as 

(y/^^o + ^{dlziYj +ffB-j^{d®d)j zi|zi j . 

For a fixed zi € R" we choose zq G 1^ in such a way that the squared bracket is equal 
to 0. Thus we obtain the estimate 

(/3Bzi|zi)-i((d®d)zi|zi)>£/3|zip, 
valid for all zi £ R". By the definition of d it holds that 

hence we obtain the identity 

PB - -{a (g) c + c (g) a) ^ I3B - -{d (g) d) + -{a (g) a + c (g) c - a (g) c ~ c Si a) 
= PB - -{d g) d) + -{a - c) (» {a - c). 

Since the matrix {a—c)g){a—c) is positive semi-definite we finally obtain the assertion. 

D 

Here is the main result on optimal L^-regularity of (|2.ip . 

Theorem 2.2. Let 1 < p < oo and assume that (H) holds true. Then (j2.ip admits a 
unique solution 

fieLp{J;H^{R"))=:Z^ 
if and only if the data is subject to the following conditions. 
(i) feLp{J;Lp{«-))^:X\ 
(ii) geLp(J;iJpl(K"))=:^^ 
(in) uo e Bpp '/^(R") =: Xp. 
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Proof. Necessity is clear by substituting the solution (w, jjl) ^ Z^ x Z^ into the equa- 
tions (123} 1 2- This yields the desired regularity for the functions /, g. The regularity 
for the initial value Uq follows from the trace theorem 

HliJ-HliR-)) n Lp(J;ff3(K")) ^ C{J-BI;^/^{W^)), 

where Bpp^^^(K") = {Hl{W^),H^{W^))i_i/p^p is the real interpolation space with 
exponent 1 — 1/p and parameter p. 

To prove sufficiency of the conditions (i)-(iii), we first apply the operator (/ — A)^^' ^ 
to both equations in (12. ip and define the new functions w = [I — A)^^/^u, rj = 

{I - A)-i/V, f ^{1- A)-i/V, ~9 = {I- A)-i/^<? and wq ^ {I - A)-i/2uo. Then 
it holds that 

/ G Lp(J;i/pi(M")), ~g E Lp{J; H^{W')), 

wo e <-'/^(M") 

and we are looking for a solution {w, rj) of the system 

wt - div(awt) = div(BV7/) + /, t > 0, a; G M", 

(2.2) rj-c-Vr]^ l3wt-Aw + g, t > 0, a; G M", 

w{0) = Wo, i = 0, a; G M", 

in the regularity class 

w G iJpi(J;7j2(M")) n Lp(J;iJp4(M")), 

77GLp(R+;ij3(M»)). 

In a next step we want to eliminate the functions g and wq. To achieve this, let w* 
be the unique solution of the problem 

/3w*-Aw* = -5, t>0,xeW\ 

w*{0) =wo, i = 0, xgM", 

with regularity 

«;* G i/pi(J;Lp(R"))nLp(J;ij2(E")), 

if and only ii g G Lp{J x M") and wq G Bpp '^(M"). Here J denotes the interval 

[0,r]. If we even have g G Lp{J;H^{W)) and wq G Bpp^/^(M") then by regularity 
theory we obtain 

w* G i/pi(J; H^{R")) n ip(J; i/p'(M")). 

The pair of functions {v, rj) = (w — w* , 77) should now solve the problem 
dtv - div{adtv) = div(SV7/) + F, t > 0, a; G R", 

(2.3) r]~c-S/r]^ /3dtv - Ad, t > 0, x e R", 

v(0) =0, t = 0, xgR", 

where i^ is defined by 

F^f + w;- div{aw;) G Lp{J; H^ (R")). 

In order to solve p.3p we take the Laplace transform in the time variable and the 
Fourier transform in the spatial variable to obtain 
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This system of algebraic equations can be written in matrix form 



A(l-z(a|0) (Bm) 
Hl3X + m (l-^(clO) 



where A € S^, > tt/2 and ^ G M" such that |A| + |^| 7^ 0. Hence the unique solution 
to these equations is given by 



1 



m(A,0 



(l-z(clC)) -{Bm 

provided 

m(A,0 :=detAf(A,0 7^0. 
To see this we consider the function m(A,^) := ?7i(A,^)/A given by 

m(A, C) - 1 - (a|0(c|0 + /5(i?eiO - *(a + c|0 + /3(i3CIOI^I '/A = ^i(0 + ^2(A, 0, 
where Z2(A,^) := /3(i3^|^)|^p/A. Let cfij — argz^; then a short computation shows 
that 

kl+Z2|>C(0i,02)(|^l| + |^2|), 

621 < TT. Here 



provided that \(f>i 



From Proposition 12.11 and the Cauchy-Schwarz inequality we obtain 



C(0i, 02) := ^ min{l, (1 + cos(0i - (/-z))'/"}. 



(a + c|e) 



<C|a + c| 



lei 



l + l^l^ 



< C\a + c\ < 00, 



hence \4>i\ < a < tt/2 for all ^ G M". Since |02| = I ai'g'A| < </> we have 

101-021 < Cr + < TT, 

provided > 7r/2 is sufficiently close to 7r/2 and this in turn yields together with 
Proposition 12.11 

|m(A,OI = \zi + Z2\ > C(|zi| + |Z2|) > C(l + ICI^ + leiVlAI) 
or equivalently 

(2.4) |m(A,OI>C(|A|(l + |en + |en. 

Observe that the converse is also true, i.e. there is a constant C > such that 

MA,c)i<c(iAi(i+ien + ien. 

In particular it holds that m{X, ^) = if and only if |A| + |^| = 0. 

Next, let wo,ui G aH^{J;Hf(W"')) n Lp{J; H^{W"-)) be the unique solutions of 

dtil - A)wo + A^uo = F-c-VF, t > 0, x e W\ 
fo(0) - 0, 
and 

dtil - A)vi + A^wi = (/ - A)i/2i;^, i > 0, x e M", 
«i(0)=0. 
The existence of vq and ui may be seen by the Dore-Venni-Theorem. It follows that 
dt{I - A)v + A\ = 5(at(/ - A) + A^)vo, 
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and 

{I-A)^/^7j = S{I-A)i(3dt-A)v„ 
where the linear operator S is defined by its Fourier-Laplace symbol 

^(^,ii — 77— p: ■ 

m(A,0 

Note that the assertion of Theorem 12.21 follows if we can show that S" is a bounded 

operator from Lp{J;Lp{W^j) to Lp(J; Lp(R")). This will be a consequence of the 

classical Mikhlin multiplier theorem and the Kalton-Weis Theorem [^, Theorem 4.5]. 

It is not difficult to show that the symbol ^(A,^) satisfies the Mikhlin condition 

(M) max|„|<[„/2]+iSupjgj5" ICl'"'|9|'5'(A,C)| < oo, 

where a G Nq is a multiindex and [s] denotes the largest integer not exceeding s G R. 
The classical Mikhlin multiplier theorem then implies that S* is a Fourier multiplier 
in Lp(]R";C) w.r.t. the variable S, and this yields a holomorphic uniformly bounded 
family {S{X)}xe^^ C 6(Lp(]R"; C)), (j) > 7r/2. By [3 Theorem 3.2] this family is also 
7?.-bounded in Lp{J; Lp(IR"; C)) (for the notion of 7?.-boundedness we refer the reader 
to [5]). Finally, since the operator dt admits a bounded ■H°°-calculus with angle 7r/2 
we obtain from [9, Theorem 4.5] the desired property of the operator S. For the 
functions u = [I — A)^/^i(; and /i = (/ — A)-'^/^?^, this yields 

u e Hl{J- i/pi(M")) n Lp(J; i/3(M")), 

as well as 

and the proof is complete. 

D 

For later purpose we need a perturbation result. To be precise we consider coefficients 
a, c and B with a small deviation from constant ones, i.e. 

a{x) = ao + ai{x) , c{x) — cq + ci{x) , B{x) ^ Bq + Bi{x), 

with ai,ci e M^i(R";]R"), Bi e VKi(R";R"^") and 

Jaijoo + Jcijoo + \Bi\oc < w. 

Furthermore we assume that divai(j;) = divci(a::) = for a.e. x S R" and that the 

quadruple (/3,ao,Co,i3o) satisfies (H). Observe that if w > is sufficiently smah, then 

(/3, a{x),c{x),B{x)) satisfy (H) as well for all x £ ^, with a possibly smaller constant 

e>0. 

We have the following result. 

Corollary 2.3. Under the above assumptions on the coefficients the statement of 
Theorem \2.2\ remains true, if uj > is sufficiently small. 

Proof. By a shift of the function u we may assume that uq = g = 0. For the time 
being we consider an interval Jg = [0,S], with a suitable small 6 > 0, to be chosen 
later. The corresponding function spaces are denoted by Xj and Zj . Moreover 

oZ} := {u e Z] : u\t^o = 0}. 

Assume that we already know a solution (w, ^) G qZ^ x Zg oi (|2.1I) . Thanks to 
Theorem 12.21 we have a solution operator S G 13{X^ x X| x Xp] qZ] x Z|) for the 
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constant coefRcient case (/3, ao,co,-Bo)- With the help of S we write the solution in 
the following way. 



_ - 




"/ 


r -\ 


u 


==5 





+ ST 


u 


.^. 









M 









where 



r 



(iw{ai{x)dtu) + div(_Bi(x)V/^) 

ci{x) • V/i 





we obtain the estimate 



(2.5) 



From the boundedness of S and since divai(a;) =0, .t G 
\{u,^^)\zlxz^, < C(\f\xl + |T(u,m)Ix|) 

for some constant C > 0. The problem is that the term |V/i| does not become small in 
Lp{Js; Lp(R")), since the function /i has no regularity w.r.t. the variable t. However, 
we have the following result. 

Proposition 2.4. Let (u,ij,) ^ Zg x Zg be a solution of (|2.ip with g = uq = 0. 
Assume furthermore that the (variable) coefficients satisfy the above assumptions. 
Then there exists a constant C > 0, independent of Jg, such that the estimate 

(2-6) |/i|i^(j,.^i(R„)) < C{\f\xi + \9\xl + HL^iJs,H^iR"))) 

is valid. 

Proof. The proof follows the lines of the proof of Proposition l3.3l D 

Owing to (|2.5I) and Proposition 12.41 we obtain the estimate 

(2.7) \iu,tJ-)\zlxZl < C(|/|xi + w|(w,M)lzixZf + |u|Lp(,/,;H|(R")))- 

The mixed derivative theorem and Sobolev embedding yield 

oH^iJs; -ffp(K")) n Lp( J,; i?3(K«)) ^ ^H^^Js; H^{M.")) -> L2p{Js; i?p(R")), 

hence by Holder's inequality we obtain \u\L^(^J^.H^^^,^)•J < S^^'^pC\u\zi and the con- 
stant C > does not depend on 5 > 0, since u|(=o — 0. Choosing first uj > 0, 
then (5 > small enough and shifting back the function u, we obtain from (|2.7p the 
estimate 

l(u,Ai)lzixz| < C'd/lx] + Ifflxf + \uo\xp), 
where C > is some constant. The latter estimates show that the operator L E 
B{Zl X Zj-X] X Xj X Xp), defined by 

dtu — div(a9tw) — div(_BV/i) 

li ~ c ■ V/i — jSdtU + Aw 

u\t=o 

is injective and has closed range, hence L is a semi-Fredholm operator. Replacing the 
coefficients (/3, a, c, B) by 

(/3, Or, Cr, Br) := (1 - T){f3, ao, Co, Bq) + t{/3, a,c,B), re [0, 1] 

we may conclude from the considerations above that for each r € [0,1] the cor- 
responding operator L-r is semi-Fredholm as well. The continuity of the Fredholm 
index yields that the index oi Li = L is 0, since Lq is an isomorphism, by Theorem 



L 
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12.21 A successive application of the above procedure yields the claim for the time 
interval J = [0,T]. The proof is complete. 

n 

3. The Linear Cahn-Hilliard-Gurtin Problem in M" 

In order to treat the case of a half space, we consider first constant coefficients which 
are subject to the following assumptions: B — bl and (a|e„) = (c|e„) = 0, where 
e„ := [0, ... ,0, —1]^ is the outer unit normal at 9M". Furthermore we assume that 
(/?, a, c, B) satisfy (H), whence it holds that 6 > e > 0. Moreover the boundary 
conditions on a and c yield that the last components of a and c are identically zero. 
We are interested to solve the following system in M". 

dtu ~ dW{adtu) = bAfi + f, t > 0, {x' , y) e M'^, 

H-c-Vfi^ (3dtu ~ Au + g, t > 0, {x' , y) e K'^, 

(3.1) dyfi = hi, t> 0, x' e M""\ 2/ = 0, 

dyu = h2, t> 0, x' £ M"~\ y = 0, 

u{0)^UQ, t = 0, {x',y)eRl. 

Note that the conormal boundary condition (W/i|e„) = hi is equivalent to ~bdy^ — 
hi, where 6 > is constant. Hence it suffices to consider the boundary condition 
dyfj, = hi with some scaled function hi. Concerning optimal Lp-regularity of p.l[) we 
have the following result. 



Theorem 3.1. Let 1 < p < cx), p ^ 3/2 and assume that (H) holds true. Then ()3.ip 
admits a unique solution 

ueH;,{J;H;,{Rl))nLp{J;H'p{Rl)) =: Z\ 

if and only if the data is subject to the following conditions. 

(i) f£Lp{J-Lp[Wl))=:X\ 
(ii) g£Lj,{J-Hl{Rl))=:X\ 

(iii) hi e Lp(J;T<-^/P(M"-i)) =: Y\ 

(iv) h2 e Wl~^'^^{J-Lp{W^-^)) n Lp(J; WV'-i/^(K"-i)) =: Y\ 

(v) UQ e Blp^'^iWl) =:Xp. 

(vi) dyUo = /i2|t=o J/P > 3/2. 

Proof. The necessity part follows from the equations and trace theory, cf. [6]. Con- 
cerning sufficiency, we first reduce (|3.ip to the case hi = h2 = uq = 0. For this 
purpose we solve the elliptic problem 



(3.2) 



(/ - A.,,)r, - d^r, = 0, x'e W'-\ y > 0, 
dyT] = hi, a;'eR""\ y ^ 0. 



Define L := {I - A^,)^/'^ in Lp(R"-i), with D{L) = H^{W^-^) and let L denote the 
natural extension of L to Lpjoc{R+;Lp{W-^)), that is D{L) = Lpjoc{R+;Hp{R''^^)) 
and Lu — Lu for each u G D{L). Then the unique solution 77 of (|3.2p is given by 
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Since hi g ip(J; Wp"^^^(R"-i)) = Dl{1 - l/p,p), we have e-^yhi e D{L) and 
therefore ry G Lp(J;_ffp(R" )), with 9y?7|y=o ~ hi. In order to remove /i2 and mq, we 
solve the initial boundary value problem 

iSdtv - A^,v -d^v^O, t> 0, x' e R""\ y > 0, 

(3.3) dyv = /i2, i > 0, x' e K""\ y = 0, 

w(0)==uo, i = 0, x' eR""\ y > 0. 

To this end we extend uq G ^pp^ (R" ) to a function uq & Bpp '''(R") and solve the 
heat equation 

I3dtv - Aw = 0, t>0, xeW\ v{0) ^uo, t = 0, xe R", 

in Lp(J;iJ^(R")). This yields a solution 

If wi := P-D denotes the restriction of v to the half space R" , the function V2 := v — vi 
should solve the initial boundary value problem 

Pdtv2 - A^,V2 - dlv2 = 0, t>0, x' £ R"-\ y > 0, 

(3.4) dyV2 = h2, t>0, x' e W'-\ y = 0, 

^;2(0) = 0, t = 0, x' £W"-\ y>0, 
where /12 := ^2 — dyVi\y^o. Set U3 = (/ — A2;/)^"i;2. Then V3 is a solution of 
/3dtV3 - A,,V3 - d^yV3 = 0, t>0, x' e R"-\ y > 0, 

(3.5) dyV3 = h3, t>0, x' e R"-\ y = 0, 

1,3(0) = 0, t = 0, a:'eR""\ y >0. 

with /13 = (/- A,0'/'^2 G oW'p^'"'^'''(^;ip(K""')) n Lp(J;iyp"'/P(R"-i)). We 
define L — {/3dt — Ax'Y^"^ with natural domain 

DiL) = oi?pi/'(J;Lp(R"-i))nLp(J;iJpi(Rr')). 
Then, the unique solution 113 of (|3.5p is given by 

i;3(y) = -L-^e-^yhs, 
and /13 G -DlII ^ ^/p,p)- This yields 

t;3 e oi^p(J;ip(K+))nLp(J;ij2(R^)). 
On the other hand, if we consider the function 114 :— dyV2 as the solution of 
I3dtvi - A^.'y4 - d^Vi = 0, t > 0, a;' e R""\ y > 0, 

(3.6) 1-4 = /i2, t > 0, x' e R"~\ y = 0, 

1,4(0) = 0, t = 0, x' e R"-\ y > 0, 

we obtain z;4(y) — e^-^^ft.2 and ft,2 G -Dl(2 — l/p,p). This yields 

z;4 e oH^{J;Lp{Rl))nLp{J;H^p{Rl)). 
From the regularity of V3 and W4 we may conclude that 

V2 e oi^p(J;i?p(K:^))nLp(J;H3(K»)). 



12 MATHIAS WILKE 

Now the functions ui := u — v and /ii :— fi ~ rj, with v = vi + V2, should solve the 
system 

dtui - div(a9tui) = 5A/ii + /i, t> 0, a;' e M"~\ y > 0, 

fii-c- V/ii = /3(9tMi - Aui +gi, t> 0, x' £ R"~\ y > 0, 

(3.7) dy^ii^O, t>0, x' eM."'-\ y = 0, 



dyUi =0, i > 0, x' e M""\ 2/ = 0, 
Mi(0) =0, i = 0, a;' e M""\ y > 0, 



with some modified data /i G X^ and gi G X^. In a next step we extend the functions 
/i and gi w.r.t. the spatial variable to M" by even reflection, i.e. we set 

. /, , ^ \fiit,x',y), iiy>0 u ' ^ }9iit,x\y), ii y > 

f2{t,x ,y) = i , and g2{t,x,y) = i a -f / n ' 

[fi[t,x',-y), iiy<0 [gi{t,x' ,-y), li y < 

Thanks to Theorem 12.21 we can solve the full space problem 

dtU2 - div{adtU2) = 6A^2 + /2, t>0, xe M", 
(3.8) fi2-c-Vfi2--l3dtU2-Au2 + g2, t > 0, x e M", 

U2(0) = 0, i = 0, xeM", 
since /2 G Lp(J;Lp(M")) and 52 G Lp{J;Hp(M."-)). This yields a unique solution 
M2 G ffp(J; i?p(K")) n Lp( J; iJp3(R")) and ^i2 e Lp(J; i?p'(M")), 

by Theorem 12. 21 At this point we emphasize that the equations p.7p ^ 2 ^^^ invariant 
w.r.t. even reflection on the hyper surface E"~^ x {0} in the normal variable y, due 
to the structure of the coefficients. This in turn implies that the solution {u2,fJ-2) is 
symmetric, w.r.t the variable y and this yields necessarily, i9j;M2|i;=o — dy^2\y=o — 0. 
Denoting by P the restriction of the solution (w2, M2) to the half space R" , it follows 
that {ui, jii) = P{u2, ^2) is the unique solution of (|3.7p and therefore u = v + ui and 
/i = ?7 + /ii is the unique solution of (|3.ip . The proof is complete. 

D 

For later purposes we will need the following perturbation result. Let Bq = bol, 
a{x) — ao + ai{x), c{x) — cq + ci{x), B{x) = Bq + Bi{x), D{x) — I + Di{x) 
with ai,ci e T4^i(R!;:;R"), Bi e M^^ (Ml: ; M" ^ " ) , Di e 14^^ (R'| ; M" '^ " ) and 

|ai|oo + |ci|oo + \Bl\oo + \Di\ao < W, 

for some uj > 0. Let furthermore divai(a::) = divci(a::) ~ for a.e. x S R" and 

{ao\iy{x)) == {ai{x)\iy{x)) = {co\iy{x)) = {ci{x)\iy{x)) = 0. 

If the constant coefficients (/3, ao, cq, So) satisfy Hypothesis (H) we have the following 
result. 

Corollary 3.2. Let 1 < p < 00, p ^ 3/2, ^ > and suppose that the data satisfies 
the conditions (i)-(v) of Theorem, \3.1\ and (DVMo|e„) = /i2|t=o J/p > 3/2. Under the 
above assumptions on the coefficients {a,c,B,D), there exists a unique solution 

u e Hl{j- Hl{RX)) n Lp{J- HliWl)), 
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of the system 

dtu - AiY{adtu) = div(BVAi) + /, t > 0, {x' , y) e M" , 
^ - {c\V^Ji) = Pdtu ~ A\v{DVu) + g, i > 0, (x', y) e 
(3.9) (-BV/i|e„) = hi, t> 0, x' e W'~\ y = 0, 

(L'Vule,^) = /i2, t > 0, a;' G M""\ j/ = 0, 
u(0) = uo, i = 0, {x\y)eRl, 
provided uj > is sufficiently small. 



^+' 



Proof. First of all, we reduce p.9p to the case wq = as follows. Extend the initial 
data uo € -Bpp (M") to some uq € Bpp (R") and solve the heat equation 

dtv - Aw = 0, i > 0, X e M", 
v{0) =uo, xeW\ 

to obtain a unique solution 

V e i/pi( J; i/pi(M")) n Lp{J; H^{R")) = ^\ 

for some interval J — [0, T]. If (w, /i) G Z^ x Z^ is a solution of p.9p . then the shifted 
function (u ~ v, fi) E qZ^ x Z^ solves p.9p with uq = and some modified functions 
/ G X-'^, g G X^ and /12 G o^'^- Observe that /,g and /12 depend only on f,g,h2 and 
the fixed function v G Z^ from above. In the sequel we will not rename the functions 
u,f,g and /i2- 
By the structure of the coefficients and by trace theory we obtain the estimate 

\{u,lJ')\zlxZl 

with a constant C > which does not depend on 5 > since u\t=Q = 0. The 
derivation of this estimate follows the lines of the proof of Corollary 12.31 The term 
|u|i {js;H^{R")) is of lower order and may be estimated by 



\u\laJs-m^AmV) < ^^'^"CH^i < S'/'^\{u,^l)\ 



zlxzi, 



hence this term may be compensated by the left side of the latter estimate if 5 > 
is small enough. If in addition w > is sufficiently small, the same is true for 
uj\{u, ij.)\zixz^- To estimate the term |V/i| in Lp{Js;Lp{M.^)), we use the following 
proposition whose proof is given in the Appendix. 

Proposition 3.3. Let (w, /i) G Zg x Zg be a solution of (13. 9p with uq — 0. Then 
there exists a constant C > 0, independent of Jg, such that the estimate 

(3-10) lMlLp(j,;_f/i(R!^)) < C'd/lx] + \9\xl + l^ily/ + \u\L^(JyM}imi))) 

is valid. 

Now the claim follows by applying a similar homotopy argument as in the proof of 
CoroUaryHSl 

D 



14 mathias wilke 

4. Bounded domains, Localization 

Let n C M" be a bounded domain with boundary dfl £ C^. In this section we solve 
the system 

dfU — div{adtu) = div(6V^) + /, t > 0, x £ Q, 

^ — c • V^ = Pdtu — Au + g, t > 0, X £ Q, 

(4.1) IN 11- v^ hi, t>0,xedn, 

dyU — h2, t > 0, a; e dfl, 
u(0) ~ uq, t = 0, a; e 51, 

with coefficients a,c G [(7^(57)]" and b E C^{il). We furthermore assume that 
div a(a::) — divc(a;) =0, x E fl, {a{x)\v{x)) = {c{x)\v{x)) =0, x E dil and {l3,a,c,b) 
satisfy (H). Before we start with the localization procedure we prove two lemmata, 
which are interesting for their own. 

Lemma 4.1. Let I < p < oo, p ^ 3/2, J — [0,T] and ft C K" be a bounded domain 
with dfl £ C^ . Then for each (3 > the initial-boundary value problem 

(idtu -Au^f, te J, X efl, 

(4.2) d^u = 5, t e J, X e dil, 

u{0) — Wo, i = 0, X e i7, 

admits a unique solution 

u e H'piJ; H'pin)) n ip(J; H^in)), 
if and only if the data are subject to the following conditions. 

(i) feLpiJ;H^in)), 

(ii) g e w^-'/^^iJ- Lj^idn)) n ip( J; w^-'^^idn)), 

(in) uo G B^p-^^P{n), 

(iv) d^uo = g\t=o, provided p > 3/2. 

Proof. The 'only if part follows from the equations and well known result in trace 
theory. Indeed, given a solution 

u e HliJ- Hlin)) n Lj,{j- Hlin)), 

of (|4.2p it follows directly that / e Lp{J; Hp{il)). Furthermore it holds that 
H'^iJ; H'^in)) n ip(J; H^in)) ^ C{J; (ilpi(fl); i?3(o))^_^/^^^) ^ c{J; i?pV'^^("))> 

by trace- and interpolation theory. Hence u(0) e Bpp (il). Finally observe that 

Vu e H^{J; Lpin)) n LpiJ; H^i^)). 
Taking the trace of Vu on dft yields 

Vulan e w^~^/^p{j; Lp{dn)) n Lp{j; w^'^/p{dn)), 

the required regularity for g. Finally, since 

w^-'/^'^iJ; Lpidn)) n Lp{j; w^-'/p{dn)) ^ c{j- Blf'^idn)), 

it follows that duu{{)) = g\t=o in case p > 3/2. To prove sufficiency of the conditions 
(i)-(iv), note that by the results of Sections [2] & [3] the unique solution of the cor- 
responding full space and half space problem to (|4.2p possess the desired regularity. 
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Then the claim for a bounded domain il C M" with dfl G C^ follows from localization, 
change of coordinates and perturbation theory, cf. [5] . 

D 



The second lemma provides maximal regularity of (j4.ll) in case a = c = and b = 1, 
the so-called viscous Cahn-Hilliard equation in its linear form. 

Lemma 4.2. Let 1 < p < oo, p ^ 3/2, J == [0,r] and 51 C M" he a bounded domain 
with dfl G C^ ■ Then for each j3 > the system 

dfU — Afi = /, t £ J, X £ Q, 
fi — (3dtu + Au = g, t £ J, X £ Q, 
(4.3) d^^ = hi, t e J, X e dfl, 

d^u — h2, t e J, X £ dfl, 
u(0) — uq, t — 0, X e i7, 
admits a unique solution 

if and only if the data are subject to the following conditions. 
(i) fGLp{J;Lj,{n)), 
(ii) geLpiJ;H^pin)), 

(iii) hi(ELp{J;Wl-'^''idn)), 

(iv) h2 e w^p"'/'P( J; Lpidn)) n Lp( J; w^-'^^idn)), 
(v) uo e B^p^/^in), 

(vi) d^uo = h2\t=o, provided p > 3/2. 
Proof. By Lemma l4. 1 1 there exists a unique solution 

V e Hl{j-Hl{n)) n L.p[j-Hl{n)). 

of the problem 

PdtV ~ Av ^ -g, t e J, X eil, 

d„v = h2, t G J, X G 9r2, 

v{0) = Mo, t = 0, X £ fl. 

Hence, w.l.o.g. we may assume g — h2 — ua = in (|4.2p . with / being replaced 
by some modified function / G Lp{J; Lp{Q)), which depends at most on the fixed 
functions / and v. 

Now we want to reduce (14. 3p to a single equation for u. Suppose that we already 
know a solution of (|4.3p . Inserting (|4.3I) -^ into (|4.3L yields the elliptic problem 

^ " /3A^ ^ l3f - Au, te J, X en, 

9,y/i — hi, t e J, X £ dfl, 

for the function /i. It is well-known that for each /3 > the latter problem admits 
a unique solution fi e Lp{J; Hp{Qj), provided {/3/ — Au) £ Lp{J;Lp{ilj) and hi £ 

Lp{J;Wp~' (dfl)). Denoting by S the corresponding solution operator, we may 
write 

^"-^ hi '^ ■ 
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Inserting this expression into (|4.3P o we obtain the problem 

Pdtu — Au = h — Su, t & J, a; e ri, 
(4.4) d^u = Q, t£ J, X e on, 

u(o) = 0, t^o, X en, 

where h :~ S{/3f,hi) and Su :~ S{Au,0). Since 5 is a bounded Hnear operator 
from Lp{J; Lp{Q)) x Lp{J; Wp^^^'^idn)) to Lp{J; H^in)) it follows that S is bounded 
and linear from Lp{J;Hp{rt)) to Lp{J]Hp{Vl)). Thanks to Lemma [4.11 there exists a 
solution operator T of (|4.2p which is a linear and bounded mapping from 



Lp{j-Hl{n)) X oW^-^'^p{j-Lp{dn))cMp{j-w^-^'p{dn)) x s^V'^^l^) 

to oHl{J; H^{n)) n Lp{J; H^{n)). With the help of T we may write 



r 



'h 




Su 





-T 










_ 0_ 



We estimate 



\r{Su,0,0)\z < C|5u|Lp(,7.ffi(o)) < C|u|Lp(j://2(n)) 

by Holder's inequality. Here the constant C > does not depend on T > 0, since the 
time traces at t = are zero. A Neumann series argument yields a unique solution 
u E Z oi (|4.4p on a (possibly) small time interval J = [0,T]. Since (|4.4[) is linear and 
invariant with respect to time shifts, the solution exists global in time. 

n 

The main result of this section reads as follows. 

Theorem 4.3. Let 1 < p < oo, p ^ 3/2, J = [0,T]. Suppose furthermore that 
a,c G [C^(f2)]", b G C^(Cl). Then (|4.ip admits a unique solution 

u e Hl{j-Hl{n))cMp{j-Hl{n)) = Z\ ^^e LpiJ;H^pin)) = Z2, 

if and only if the data are subject to the following conditions. 

(i) feLpiJ;Lpin))^X\ 
(ii) geLp{J;Hl{n))^X^, 

(iii) h,eLp{J;W;,-'^''{T))=Y\ 

(iv) h2 e w^-'/^PiJ;Lpir)) n Lp{j;W^-'^Pir)) - y^, 
(v) woe5pp'^''(f^) = ^p, 

(vi) d^uo = /i2|t=o i/p > 3/2. 

Proof. By Lemma [4.21 we may first reduce (|4.1|) to the case /ii = /12 = uq = and 
some modified functions f,g in the right regularity classes. We cover Q by finitely 
many open sets Uk, k = I, ..., N, which are subject to the following conditions. 

(i) C/fe n F = and Uk = Br, (xk) for ah fc = 1, ..., A^i; 
(ii) C/fc n r 7^ for fc = A^i + 1, ..., N. 
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We choose next a partition of unity {fk}k=i such that J2k=i 'fk{x) = 1 on fi, < 
fk{x) < 1 and supp (pk C Uj. Note that {u, /i) is a sohition of (I4.ip if and only if 

StWfc - diviadtUk) = div(6VAife) + fk + Fkiu^fi), t e [0,6], xennUk, I < k < N, 

^ik~c-Vfik = (3dtUk - Auk + gk + Gk{u, fi), te[0,S], xennUk, l<k<N 

(4.5) foV^ife • i^ = (6V(^fe • z^)Ai, te[0,S], xernUk, Ni + l<k<N 

duUk = ua^(/3fe, < e [0, 5], xeVC^Uk, Ni + l<k<N 

uk{o) = 0, i = o, xe rinc/fc. 

Here we have set Uk = wpk, t^k = fJ-^k, fk = f^k, gk = f^k- The terms Fk{u,^i) and 
Gk{u, /i) are defined by 

Fk{u, fi) = -(a • Vipk)dtu - (V6 • Vipk)^^ - 2bVipk ■ V^ - b^Aipk, 
and 

Gk{u, fj.) = -{c ■ Vipk)lJ' + 2VuS/tpk + uAifik- 
In case A; = 1, ..., A^i we have no boundary conditions, i.e. we only have to consider 
the first two equations in (|4.5p . In order to treat these local problems with the help 
of Corollary 12.31 we extend the coefficients from Br^. [xk ) to R" in such a way that 
diva(x) — divc(a;) = 0, x € R", holds for the extended coefficients a and c. Note 
that w.l.o.g. we may assume Xk — 0. This follows by a translation in M". 
We use the following extension a of a (or c of c). 

' a{x), 



(4.6) 



T^k 



a'=(x) = 



X(.BrM. 

X G 



\BrM. 



where r — \x\, ^ — x/\x\ and S,j, Uj denote the components of £, and a, respectively. 
The task is to compute the scalar valued function R{r,^). Since diva(a;) =0, x G fl, 

the divergence of a ( -^x) and [^ a [ ^^j ) ■? ni^Y be computed to the result 



div 



2SE^^^^^^«* 



and 



div 



3-))e 



(n-1) 



«j 



$]]-5y:e.6ft:a, 



The divergence of the last term R{r, ^)f is given by 

n- 1 



div[R{r,C)C]=drRir,0 



r-2 /L^ 



R{r,0. 



Finally, this yields that diva'^(x) = if and only if the function R — R{r,^) solves 
the ordinary differential equation 



drRir,0 



(n-1) 



Rir,0 



(n-1) 



r 



r >rk. 



In order to achieve a'',c'' S W^ ( 



initial condition R(rk,^) 
given by 



K"), we require a}'{rkC) = o.{fkS,)- This yields the 
2(a(?'fcOIC)' hence the function R ~ R{r,£,) is explicitly 



„n-l 



Rir,0-:^R{rk,0 + 



2(n- 1) 



n-2 



S^ 



£.] ds, r >rk. 
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Since 



2^(«(^fce)IO = 2(a(rfcOIO 



,(^-1) 



s''-^(a{rkO\0 ds, 



we may write 



^ -2U 



2{n- 1) 



?« 



- a{rkO 



a^ds, 



in case \x\ > r^- Owing to this identity and the assumption a,c £ C^{Vl), it is evident 
that there holds 



a'=(x) 



+ |a'=(x)-c(0)| <uj, 



(4.7) 



b''{x) 



for ah X € M", where uj > can be made as smah as we wish, by decreasing the radius 

rfc of the charts Uk, k e {l,...,Ni}. 

For the coefficient function b we use the reflection method from [5], i.e. we set 

'b{x), xeB^jO), 

It may be readily checked that b'' G W^{R") and that 

\biO)-b''{x)\<ij, xeR", 

with the same oj > as above. Hence for each chart Uk, fc e {1, ..., A^i} we have coef- 
ficients which fit into the setting of Corollarv 12.31 Therefore we obtain corresponding 
solution operators S^ of (|4.5p such that 



(4.8) 



Si 



fk + Fk{u,n) 
9k + Gk{u,n)_ 

for each /c S {1, . . . , A^i}. 

For the remaining charts Uk, k € {A^i + l, . . . , N} we obtain problems in perturbed half 
spaces with inhomogeneous Neumann boundary conditions. For the further analysis 
we have to understand how to treat (|4.ip in such a setting. To this end we fix a point 
xq € on and a chart U{xo) n Oil ^ 0. After a composition of a translation and a 
rotation in R", we may assume that xq = and I'ixo) = [0, . . . , 0, —1] = e^. Consider 
a graph p € C'^(IR"~^), having compact support, such that 

{{x',Xn) e U{xo) C M" -.Xn = p{x')} = 8^ H U (xo) ■ 

Note that by decreasing the size of the charts we may assume that jV^'ploo is as small 
as we like, since S/x'p{0) — 0. 

For the time being, we only know that diva(a;) = divc(a;) — for all x e U{xo) D fl. 
So we have to extend the coefficients a and c in a suitable way. To this end we first 
transform the crooked boundary U{xo) dft to a straight line in R"~^ x {0}. This 
will be done with the help of a suitable transformation. Let u(x',x„) = v{g{x)) = 
v{x',Xn - p{x')) and p{x) = r]{g{x)) = r](x',Xn - p{x')) and ^^^(a^o) = g{U{xo))- 
Then the differential operators a ■ Vu and c • V/i transform as follows. 

a{x)-Vu{x) ^ a{x)-iDg{x)'^Vv{g{x))) = iDg{x)a{x))-Vv{g{x)) ^ aig{x))-Vv{g{x)), 

and 

c{x)-W^i{x)=c{x)-{Dg'^{x)Vr,{g{x))) = {Dg{x)c{x))-VT^{g{x)) = c{g{x))-VT^{g{x)), 
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with a{x) :— Dg{x)a{g^^[x)) and c{x) — Dg{x)c{g^^(x)). The transformed Laplace 
operator reads 

Au = d\v{DgDg'^Vv). 
Similarly we obtain 

div(6V/i) = div(i?V7/), 

where B{x) := b{g^^{x))Dg{x)Dg'^ (x), x e g{U{xQ)) nWl- Here the matrix Dg is 
given by 

In^l 0" 

_-V^,,p{x'y I 

where /„-i is the identity matrix in Mf"^^^^'"^^^ Observe that the normal v at 
U{xo) n dn is given by 

iy{x',p{x')) 



Dg{x) 



x' e 



n/T+F^ 



-1 



Therefore it holds that ^JlTW^^M^'WiDg'^y^v = [0, ... ,0, -1]"^ = e„, hence the 
transformed boundary conditions are JSVry • e„ = y/1 + \\7p{x')\'^Q^^hi and 



DgDg'^Wv ■ e„ = y^l + \Vp{x')\'^e-^h2. 

Here O^^ is defined by {Q~^u){x) := ■(/(^"^(x)), x e K!;:. 

By construction, the transformed coefficients satisfy diva(x) = divc(x) = for all 

X G Broixo) n R^ and (a(a;)|e„) = (c(a;)|e„) = for ah x G Br-oixo) n SM!;.. Now we 

are in a position to use the extension (14.61) in order to extend a and c to the whole 

of K", such that the divergence condition diva(a;) = divc(a;) = is preserved for 

X G R" . It is furthermore clear by the structure of ()4.6p that (a(a;)|e„) = (c(a;)|e„) = 

holds for all x G dM.^ — M"^^ x {0}. The coefficient matrix B can be extended 

to a matrix B on M" by the reflection method (|4.7p . In particular it holds that 

B(xq) = B{xo) = B{xq) = b{xo)I, by construction. 

Therefore wc have to solve the following perturbed problem in the half space M" . 

dtv - div{ddtv) = div(i?V77) + O^V, te[0,S], xe M" , 

ri-c-Vr]^ f3dtv~diviDVv)+e-^g, te [0,6], x G R" , 

(4.9) (SV77|e„) = Vl + |V.'p(x')Pe-i/ii, t G [0,<5], x' G R"-\ y = 0, 



pVi;|e„) = v/l + |V,'p(a;')Pe-i/i2, t G [0, <5], x' G R"-\ y = 0, 
?;(0) = 0, t = 0, a; G R", 

withD := DgDg'^ G I^^(M""i) and some functions (/, 5, /ii, /la) e X^xX'^xY^xY'^ 
such that ft.2|t=o — 0. From the extension method above it follows that 

\a{x) - a{xo)\ + \£{x) - c{xq)\ + \B{x) - B{xo)\ < uj, 

for all a; G R" where we can choose u > arbitrarily small, by decreasing the radius 
To > of the ball Brg{xo) — g{U{xo))- Furthermore it holds that \D{x) ~ I\ < ut, x £ 
R^, since we may choose |Vp|oo as small as we wish. An application of Corollarv l3.2l 
yields a unique solution operator S^ of (|4.9I) . hence QS^ is the corresponding solution 
operator for the chart U{x(j). At this point we want to remark that the function 
VTT]WpF is a multiplier for the spaces Wl~^'^{W-^) and Wp~^'^{W~^), since 
p G C'^(R"~^) has compact support. 
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This above computation yields solution operators &kS^ for the charts Uk, k £ {Ni 



1, . . . , N}, hence we may write 
(4.10) 



Gfe^f 



'ek\^ + Fk{u,^i))■ 

Qk\9k + Gkiu,n)) 
Ql'^{ud„ipk) 



for each k G {iVi + 1, 
{1, . . . , A^}, we obtain 



.,N}. Summing (I4.8p and (|4.10p over all charts Uk, k £ 



(4.11) 



k=l 



fk + Fk{u,^) 
9k + Gk{u,fj.) 



N 

E 



65,^ 



'&-'{fk+Fkiu,fi))' 

e-\gk + Gk{u,n)) 

e-i(SV^fc-i^)Ai 

Q^^iud^ifk) 



since {"ySfclfcLi is a partition of unity. By the boundedness of the solution operators 
we obtain the estimate 

(4.12) 

l(", Ai)lzi xZ| < M{\f\xi + \9\x] + l"|Lp(Jo;ff|(n)) + |9tM|Lp(Jo;Lp(n)) + |/^|Lp(Jo;ffi(n))), 

for some constant M > which is independent of the interval Jq — [0, 6] under 
consideration. The term \u\i^ {Jo;H^{n)) may be estimated by S^''^^C\u\zi with some 
constant C > being independent of Jq. To estimate the remaining terms we need 
the following result. 

Proposition 4.4. There exists a constant M > 0, independent of Jq, such that 

lAi|Lp(Jo;ffi(0)) + |5tw|Lp(Jo;Lp(n)) < ^d/lx] + Ifflxf + l^llly/ + I"|Lp(,/o;//|(0)))- 

Proof. The proof follows the lines of the proof of Proposition 13.31 D 



Choosing (5 > sufficiently small, we obtain from (I4.12p and Proposition 14.41 the 
estimate 



Xlx 



\iu,^J.)\zyKZ| < M{\f\xi + Iglx] + l^ily/ + I'i2|y/ + \uq\xJ, 

for a solution of (|4.ip . This shows that the bounded operator L : Zg x Zg 
Xg X ys defined by 

dtu — div{adtu) — div(i3V/i) 
l-i — {c- V/i) — I3dtu + Am 
L{u,fi)= [BV^-v) 

is injective and has closed range, i.e. it is semi Fredholm. Here y^ is defined by 
ys := {(/ii,/i2,Mo) e r/ X y/ x Xp -. d,uo = /izk^o, P > 3/2}, 

which is a closed linear subspace of the Banach space Yg x Yg x Xp. To show 
surjectivity, we apply again the Fredholm argument to the set of data 

(/3,a,,c„B,)-(l-T)(/3,0,0,/„)+r(/3,a,c,B), re [0,1]. 

The corresponding operators Lr are semi Fredholm by the above procedure and by 
Lemma 14.21 the operator Lq is bijective. The continuity of the Fredholm index thus 
yields that the index of ii = L is and therefore the operator L is bijective as well. 
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A successive application of the above arguments yields existence of a unique solution 
(m, /i) of ()4.ip on an arbitrary bounded interval [0,T]. This completes the proof of 
Theorem [131 

D 

5. Local Well-Posedness 

Let p > n + 2, / e X\ g g X'^, hj e Y^ , j = 1, 2 and Vo & Xp be given such that 
the compatibility condition d,^^o = /i2|t=o is satisfied. In this section we consider the 
quasilinear system 

dtil^ - div(a(x, V', VV')5tV') = div(6(x, ip, V'0)V^) + /, t > 0, x e ft, 

fi - c{x, ip, VTp) ■ V/i = Pdfil^ - A-0 + $'(-0) +3, t > 0, a; e ri, 

(5.1) h(x, ip, Vip)d^n = /ii, i > 0, a; e r, 

d^tp = /i2, i > 0, X e r, 

V-W = V'o, t = o, X en, 

where $ G C'^^(M). Assume that we have given vector fields a,c € C^{fl;C^^{R x 

M":M")) and a scalar valued function b G C^{n- C'^-{M. x M";]R)) such that 

(5.2) 

a{x) := a{x,ijo{x),Vxpo{x)), b{x) := b{x,i>o{x),Vi>o{x)), c{x) := c{x,4>o{x),'^4>o{x)) 

satisfy the conditions 

(5.3) diva(a;) = divc(a::) = 0, x G fJ, 



(5.4) 



mx)\iy{x)) = {dix)\iyix)) = 0, a; G dfl. 



Suppose furthermore that (/3, a, c, b) are subject to Hypothesis (H) for each a; G 51. Ob- 
serve that for p > n+2 wehaveV'o G ATp = Bpp^^'\n) ^ (7^(11), hence 5, 2 G [C^{Tl)]" 
and b G C^(fl) and therefore the coefficients, frozen at ■00, satisfy the assumptions in 
Theorem nil 

Thanks to Theorem 14.31 we may define a pair of functions {u*,v*) E Z^ x Z'^ as the 
unique solution of the linearized system 

Uj - div(au*) = div(6Vw*) + /, t>0, x en, 



(5.5) 



We set 



- 6 . Vw* = I3uf - Am* + g, t > 0, x e n, 
bVv* -v ^hi, t > 0, a; G r, 
d^u* = ft,2, t > 0, a; G r, 
M*(0) = Vo, i = 0, a:G f^. 



El = Z^{T) X ^2(r)^ pEi = {{u,v) e El : u|t=o = 0}, 
Eo = X\T) X X\T) X yi(r) X Y^{T), qEq = {(/,ff, /ii,/i2) e Eo : /i2|t=o = 0} 
and denote by | • |i and | • |o the canonical norms in Ei and Eq, respectively. We define 
a linear operator L : Ei — >■ Eq by 

dtu — (lYv{adtu) — div(6Vw) 

'} — c ■ Vw — pdtu + Au 

b'S/v ■ V 

dvU 



L(u, v) = 
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and a nonlinear function G : oEi x Ei — > qEo by 

G{(u,v),{u\v*)) 



Gi{{u,v),{u\v*)) + G2{{u,v),{u\v*)) 

G^{{u,v)Au\v*)) + Gi{{u,v),{u\v*)) 

G^{{u,v),{u*,v*)) 

0, 



where 



Gi[u,u*) = div[(a(a;,u + u*,V(u + u*)) - a)dt{u + u*)], 

G2{{u, v), {u\v*)) = div[(6(2;, u + u\V{u + u*)) - b)V{v + w*)], 

Gz{{u,v),{u*,v*)) = {c{x,u + u*.y{u + u*)) -c) • V(t' + w*), 

Gi{u,u*) = $'(u + u*), 



and 



G5((u, w), [u*,v*)) ^[b- bix, u + u\V{u + u*))]V{v + v*) ■ v. 

Considering L as an operator from qIEi to qEq, we obtain from Theorem 14.31 that L 
is a bounded isomorphism and by the open mapping theorem L is invertible with 
bounded inverse L"-'^. It is easily seen that (V', [i) := (w + m*, w + w*) is a solution of 
(|5.ip if and only if 

L(u, v) — G{{u, v), {u* ,v*)) or equivalently (u, v) — h^^G{{u, v), {u* , v*j). 

Consider a ball Mr C qEi where r e (0, 1] will be fixed later. Define a nonlinear 
operator by T{u,v) :— L~^G((u, w), (u*, w*)). To apply the contraction mapping 
principle we have to show that TTBr C B^ and that there exists a constant k < 1 such 
that the contractive inequality 



(5.6) 



|T(m,u) -T{u,v)\i < k\{u,v) - {u,v)\i 



holds for all {u,v), {u,v) € B^. The following proposition is crucial to prove the 
desired properties of the operator T- 

Proposition 5.1. Let p > n + 2, J = [0,T] and assume $ G C^~(IR). Then there 
exists a constant C > 0, independent of T and r, and functions /ij — ^J^j(T) with 
fJ-j{T) — > as r — > 0, j = 1, . . . , 5 such that for all (mi, wi), (u2, ^2) € Br the following 
statements hold. 

(i) \Gi{ui,u*)-Giiu2,u*)\xi < C'(r + Mi(T))|(wi,{)i) - {u2,V2)\i; 
(ii) \G2iiui,ii)Au*,v*)) - G2iiu2,i2),iu*,v*))\x^ < C(r + ^2 (T)) | (ui , i^i) - 

("2,w2)|i; 

(iii) \G3iiui,vi),{u*,v*)) - G3{{U2,V2)AU*,V*))\X- < C(r + /i3(r))|({ii,«i)- 

{u2,V2)\i; 
(iv) \G4:{U1,U*) ~ G4:{u2,U*)\x2 < Cp4:{T)\{ui,Vi) - iu2,V2)\i; 

(v) \G5iiui,ii)Au*,v*)) - G5iiu2,V2),iu*,v*))\yi < G{r + p5iT))\iui,vi) - 

{u2,V2)\l- 

Proof Define the ball Br(u*, w*) C Ei by means of 

Mriu*,v*) := {(u,u) e El : (u,w) = {u,i) + (m*,w*), (u^v) £ BJ. 
Let {uj,Vj) G Br{u*,v*), j e {1,2}. Observe that 

\uj - u*\oo,Xp < Co\uj - u*\zi- < r, 
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with some Co > 0, which is independent of T > 0. This yields 

\Uj\oo,Xp <Cr+ |m*|oo,Xp <Co + \u*\oo,Xp =■ R, 

since r S (0, 1]. To prove the first part, note that 

div[(a(a;, mi, Vmi) — a)dtUi] — div[(a(x, M2, VU2) — a)dtU2] 

— (a(x, Ml, Vui) — a) • "S/dtUi — {a{x, U2, VM2) — a) ■ 'S/dtU2 
+ div(a(a;, ui, Vmi) — a)dtUi — div(a(a;, U2, VM2) — a)dtU2- 

Next we have 

{a{x, ui, Vui) — a) ■ VdtUi — (a(x, U2, VU2) — o.) ■ V9fU2 
= {a{x, ui, V-ui) - a(x, U2, VM2)) • V9tMi + (a(a;, M2, VM2) - a) ■ {VdtUi - VdtU2)- 
Therefore we may estimate 

|(a(-,-ui, V-ui) - a(-,M2, VU2)) • Vdtui\xi 

< |a(-,ui, Vmi) - a(-,U2, Vu2)|oo,oo(|V5fWi - 'Vdtu*\x^ + \'Vdtu*\xi) 

< L{R)C{r + A*(r))(|Ml - M2|oc,oc + |Vmi - Vu2|oo,oo) 

<L{R)C{r + ^x{T))\ul-U2\z^, 
as weh as 

|(a(-,U2, VU2) -a) • iydtui ~ VdtU2)\x^ 

< (|a(-,Mi, Vmi) - a(-,M*, Vm*)|oo,oo + |a(-,M*, Vm*) - a|oc.,oo)(|V9tUi - VdtU2\x^ 

<L{R)Cir + fi{T))\ui~U2\zi, 

where ^{T) := max{|V9fU*|xi, \u* — ipo\oo,Xp} — ;► as T -^ since u* e Z^ is fixed 
and u*\t=o — ipo- For the remaining terms we use the identity 

(5.7) div(a(x, u, Vu)) — divj,a(a;, u, Vu) + dza{u, Vw) • Vu + dqa{u, Vu) : V^u, 

where a = a{x, z,q), q = [qi, . . . , qn]^ 

n 

dqa{u,\/u) ; V^u := >, dq.aj{u,'^u)didjU. 

Furthermore we make use of 

div(a(a;, ui, Vui) — d)dtUi — div(a(x, U2, VU2) — d)dtU2 
= (div(a(a;, ui, Vui))— div(a(a:, U2, Vu2)))9tWi+div(a(a:, U2, \7u2)~a){dtUi—dtU2)- 

Let us first estimate div(a(a;, ui, Vui)) — div(a(x, U2, VU2)) in ioo(0, T; Loo(^))- By 
(15.71) we obtain 



|div(a(-,Mi,Vui)) - div(a(-,M2,Vu2))|oo,oo 

< \dza{-,Ui,Vui) - 9za(-,M2,Vu2)|oo,oo|Vui|oo,oo + |92a(-,W2,Vu2)|oo,oo|Vui - Vu2|oo,oo 

+ \dga{-,ui,Vui) - 9ga(-,U2, Vu2)|oo,oo|V^ui|oo,oo + |9ga(-, M2, Vw2)|oo,oo|V^ui - V^W2|o 

< C{R,u*)\ui - W2I1, 

where C{R, u*) > and depends only on R and the fixed function u* £ Z^ but not 
on T and r; recall that r S (0, 1] and (wi — U2)|t=o = 0. Furthermore 

\dtui\xi < \dtui - dtu*\xi + \dtu*\xi < C{r + |5,u*U0, 
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with C > being independent of T and |9tw*|xi — >■ as T — ;► 0. In a similar way we 
obtain 

|div(a(-,-U2, Vm2)) - div(a(-,?Ao, VV'o))|oo,oo 

< |a2a(-,-U2, Vua) - 92a(-, V'o, VV'o)|oo,oo|V?/'o|oo + |92a(-, W2, Vu2)|oo,oo|Vu2 - VV'olocoo 

+ |9,a(-,M2,VU2) -9qa(-,'0O, VV'o)|oo,oc.|V^Vo|oo + |9ga(-,M2, Vu2)|oo,oo|V^U2 - V^Vo|oo,oo 

< C{R,U*)\U2 --(Aoloo.Xp 

<C{R,u*){r + \u*-Mo.,x,). 

Note that \u* — ?Ao|oo,Xp — ^ as T — > since w*|t=o = V'o- FinaUy it holds that 
\dtUi — 9fM2|xi '^ Wi ~ ■'^2|i- This proves (i). Statements (ii) and (iii) follow in a very 
similar way, while (v) follows from trace theory and (ii). To prove (iv), we use the 
condition $ € C^~ (R) to conclude 

|$'(ui) - $'(U2)U2 < T^/P{W{ui) - $'(«2)|oo,oo + |V$'(ui) - V$'(7.2)|oo,oo 

< T^'PC{R){\ui - M2|oo,oo + |$"(wi)|oo,oo|Vwi - Vu2|oo,oo + 1^2 |oo,oo | $"(ui) ~ $"("2) | 00,00) 

< T^/PCiR, U*){\U1 - U2\oo,o. + |Vmi - Vw2|oo,oo) 

<T^/PC{R,u*)\ui-U2\i, 

where C{R, u*) > does not depend on T > and r e (0, 1]. The proof is complete. 

D 



With the help of Proposition 15.11 we are able to prove the desired properties of the 
operator T defined above. We first care about the contraction mapping property. 

(5.8) 

|T(mi,Wi) - T(u2,W2)|l < |L"^||G((ui,Ui), {u*,V*)) - G{{U2,V2), {u*,v*))\o 

<|L-1|(|Gi(mi,7/*)-Gi(m2,w*)|xi 

+ |G2((mi,Ui), (u*,W*)) - G2((W2,W2), {u*,V*))\xi 

+ |G3((ui, i;i), {u*,v*)) - G3iiu2,V2), iu*,v*))\x2 

+ \G4{ui,U*) - G4{u2,U*)\x2 

+ \G5{{ui,v,),{u*,v*))-G5{{u2,V2),{u*,v*))\y^) 

< C{r + fi{T))\{ui,vi) - {u2,V2)\i. 

where ^ = ^(T) is a function with the property that fJ,{T) — ?► as T — ^ and G > 
is a constant which does not depend on T > 0. Thus, if T > and r g (0, 1] are 
sufficiently small we obtain ()5.6p . The self mapping property can be shown in a similar 
way. The above computation yields 

ir(u,z;)|i<|r(^i,«)-r(o,o)|i + |r(o,o)|i 

(5.9) <G{{r + fiiTmu,v)U + \G{{0,0),{u*,v*j)\o) 

<G{{r + ^i(T))r+\G{{0,0),{u*,v*))\o). 

Since G((0, 0), {u*,v*)) is a fixed function in Eq it follows that |G((0, 0), {u*,v*))\o -^ 
as r — >■ 0, whence TTB^ C B^, provided that T > and r G (0, 1] are small enough. The 
contraction mapping principle yields a unique fixed point (u, v) £ qKi or equivalently 
(ipjfj.) := {u + u*,v + v*) G El is the unique local solution of ()5.ip . Therefore we have 
the following result. 
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Theorem 5.2. Let p > n + 2, Jo = [0,To] and suppose that $ e C^^(R), a,c e 
CH^;C2-(RxM";R")) andbe C^^; C2-(MxR"; M)). Then there exists an interval 
J = [0,T] C Jo, suc/i t/iat (j5.ip admits a unique solution 

if the data are subject to the following conditions. 

(i) feLp{J;Lp{n))=X\ 
(ii) geLp{J;H^p{n))^X^, 

(iii) /iieip(j;W-;-^/^(r)) = yi, 

(iv) /i2 e w^p^'/'^C J; £p(r)) n ip( J; w^-'^^iT)) - y2, 

(v) Vo e s;^p"'/''(^) = ^p, 

(vi) 9^-00 = h2\t=o, 
(vii) (/3, a, 6, c) satisfy (H) for all x d Vt as well as (I5.3p and 



Remark 5.3. ^n inspection of the proof of Theorem \5.S\ shows that the assumption 
p > n + 2 can be relaxed to p > {n + 2)/3 in the semilinear case, i.e. if (a, &, c) are 
independent of ^p and Vijj. Indeed, it remains to estimate the nonlinearity $'("0) in 
Lp{0,T; Hp{Q,)). However, in the sequel we will always assume the stronger condition 
p> n + 2. 

6. Global Well-Posedness 



Let n < 3 and p > n + 2 according to Theorem 15.21 In this section we consider the 
semihnear version of (|5.ip . i.e. we assume that a — a{x), c = c{x) and B = b{x)I. 
Then, a successive apphcation of Theorem 15.21 yields a maximal interval of existence 
■/max — [0,T,„ax) for the solution {ipTpi) of (|5.ip . i.e. (|5.ip admits a unique solution 
and 

^^Hl{J■Hl{n))f^Lp{J■Hl{n)), fieLp{J;H^p{n)), 

for each interval J = [0,T] C Jmax- 

Suppose T,„ax < oo and let J — [0,T] C [OjTmax)- We start with an a priori estimate 
for the solution ■;/; G Z^ on the maximal interval of existence Jmax- To do so we 
multiply (|5.1P ;^ by /i, ()5.iP o by —9(0 and integrate by parts to obtain 

(6.1) / (^t^p^i + (BVAi|V/i) + {a\\'^i)^t^^ dx =^ nf dx + ^hi dV 

and 
(6.2) 

J (-at0A^+(c|V/i)9t0+/3|at0|2 + i^|VVP + ^<i>(^)) dx = J Mh2 dV-j dti^gdx, 
since (a|t^) = on dO,. Adding (|6.ip and (|6.2p yields the equation 

(6.3) Jt(\\'^i'\l + j^ '^(^) dA +I3\dti:\l + {a + c\dt^pVfi)2 + (SV^I Va*)2 

= fif dx + iihi dT + dtiph2 dT ~ dttpg dx. 
Jn Jr Jr Jn 

From Assumption (H) with zq — dtip and zi — V/i it follows that 

'tVia + (a + c\dti^Wfi)2 + (BVm|Vm)2 > e{\dtij\l + \^t^\l)- 
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For the first and the second integral in ()6.3|) we apply Holder's inequality as well as 
the Poincare-Wirtinger inequality to obtain 

f^lfdx<C\f\2(\'^^i\2 + \f^^dx\] and f ^lhldT<C\hl\2,^(\'^^^\2 + \ f ^^dx\ 

The integral /^ /^ d,^ can be computed in the following way. Since div c = in J7 and 
{c\iy) = on r we have 

/ (c|V/i) dx = {c\i^)fJ' dr — / /idivc dx = 0, 
Jn Jr Jn 

hence it follows from (JS.IP -^, (|5.iP n and the boundary conditions that 



II dx = P / dt4' dx + $'("0) dx + / g dx 
Jn Jn Jn 

= / $'(V') dx+ g dx + f3( f dx + hi dV 
Jn Jn \Jn Jr 

Assume in addition 

(6.4) ^s)>-^s^-co, seK, 

where Cq > and < 7y < Ai, with Ai > being the first nontrivial eigenvalue of the 
negative Neumann Laplacian and 

(6.5) |$'(s)| < (ci$(s) + C2s2 + C3)^ for all s G M, 
and some constants Ci > 0, 6 £ (0, 1). This yields 



\ lidx\< (ci$(^) + C2|V^r + ca)" dx + c(|5|i + |/ii|i,r + |/|i). 
Jn Jn 

By the last estimate, Young's inequality and the Poincare inequality it holds that 
(6.6) 

l^fifdx + l^fihidT<C{S)(^\V^\l + l^<^Wdx+\f\l + \hi\l^ + \g\l + l^+S\Vii\l 

where q :— max{2, j^} and S > may be arbitrarily small. For the term J^ dtipg dx 
in (j6.3|) we apply Young's inequality one more time to obtain 



(6.7) / dt4'gdx<d\dt^\i + Ci6)\g\i 

Jn 

Integrating (|6.3p with respect to t and choosing S > small enough, we obtain together 

with (|6.6p and (|6.7I) the estimate 

(6.8) i|VV(i)l2 + / ^m)) dx + Ci{\dt^P\l2 + |VAi|^,2) 



C2 ^f^ (^|V^(r)|2 + $(V.(r))) dr + \f\l, + \hi\l^,. + \g\l, + l) 



dtijjh'2 dT dr. 
Jr 

In order to treat the last double integral, we have to assume more regularity for the 

function h2- To be precise, we assume that 

h2eH'p{J;LpiT))nLpiJ;W^''/PiT))^CiJ;Lp{T)). 
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Due to this fact, we may integrate the last term in (|6.8|) by parts to the result 

(6.9) / / dt^h2 dT dT= f tjj{t)h2(t) dV- ( Vo/t2|t=o dT - f f M^2 dV dr, 
Jo Jr Jr Jr Jo Jr 

where we also made use of Fubini's theorem. For the first term we use Young's 
inequality, the embedding H2{^) ^^ -^2(r) and the fact that 

(6.10) / V(i) dx^ tPo dx+ f dx dT+ hi dT dr. 
JO. Ja Jo Jn Jo Jr 

This yields 



mh2{t) dV < <5|V(t)|^i(o) +C(<5)|/»2(t)llr 

< SC\V^Pit)\l + C{S) (|/i2|L.2,r + l/ki + \hi\i.i,r + |V'o|i) • 
Observe that we have ft.2|t=o = d^'4'o G Bpp (r) ^-> -^2(r) and, by trace theory, 

B^p-yp{n) ^ i?pV'/^(r) -^ UT). 

It follows that the integral /p '00^2 |t=o ^r converges. Finally, concerning the last term 
in (16. 9p we apply Young's inequality one more time to the result 

J JtPdth2 dTdT<\j |V'(t)|?,1(o) dT+]^\dMl2,T 

<C f \Wi;{T)\l dT + C(TJ,hi,dth2,^o), 
Jo 

where we used again (|6.10l) . Set 

E{u) = -\Wu\l + / $(ii) dx, u e Hl{Vi). 
2 Jn 

Then by the above estimates there exist some constants Cj > such that 
Em)) + Cli\^t^\l2 + \V^l\l2)<C2 I Ei^ir)) dT + C3iTJ,g,hi,h2,dth2,^o), 



for all t e [0,r], provided that J > is sufficiently small. With the help of (|0| 
it follows that E{u) is bounded from below for all u g H2{^), hence we may apply 
Gronwall's lemma to the result that E{ip{-)) is bounded on J,„ax = [0, T„iax). Applying 
(16.41) one more time and using the fact that | J^^ ?/'(i, x) dx\ < C it holds that 

^e ioo(Jmax;i?2(^)). 

Note that in the semilinear case the following estimate for the maximal solution {ip, fi) 
of dO holds 



(6.11) \^\z^{T) + \tJ-\z^(T) 

< C (|$'(V')U2(T) + l/Ui(T) + \9\x^{T) + I^i|y1(T) + \h2\YHT) + IV'oUp) ■ 

Here the constant C > does not depend on T € (0,rmax). Suppose that ^'{tp) 
satisfies the estimate 



.12) 



I^'WIXHT) < C{TMzHT)mL^iO,T^^^;H^mr 
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for some k G (0, 1) and m > 0, where C{T) > and sup';-g[g j.^^^) C{T) < cx). 
Substituting (p7n|) into (pUj) yields 

|V'|zmt)<m(i + |V'||i(3.)), 

where M > does not depend on T e (0, Tmax)-This in turn yields that |V'lzi(T„ax) i^ 
bounded, since n € (0, 1). Therefore '0(Tmax) € Bpp (fl) is well-defined and we may 
continue the maximal solution {ijj,fj-) beyond the point Tmax, which is a contradiction 
to the maximality of Tmax- 

It remains to show the validity of (j6.12p . We start with the term W^'{ip) = ^"{ipy^ip 
in Lpiny. It holds that 

|*"WV^|p < |<i>"(V')l3p/2|VV|3p, 

by Holder's inequality. Assume that there exists a constant C > such that 

(6.13) |$"(.)|<C(l + |.r), 

for all s G ffi and some a > 1, where a < 4 in case n = 3. Then we have 

Applying the Gagliardo-Nirenberg interpolation inequality we obtain 

H\zc.p,2<C\i^\Hl(n)Wr. 

provided 

n 2n f n n\ , ^ 

- - :— = a 3 - - + - , a G [0, 1 . 
q Zap V p qj 

On the other side we obtain 
provided 

l- — + - = b{3-- + -), &e [1/3,1]. 

3p q \ P qJ 

Chose q in such a way, that -ff2(^) ^^ Lq{il), i.e. n/q > n/2 — 1. Thus q may be 
arbitrarily large if n e {1, 2} and g < 6 in case n = 3. If n = 3, let 

/^ ^ A\ an . ( 3ap\ 

(6.14) _<(^<mm|6,^|, 

while in case n = 1 , 2 we require 

an Sap 

(6.15) T<^<^- 

This is possible, since n < 3p for n < 3 and an/2 < 6 if n = 3, since in this case we 
assume a < A. Now it follows that 

(6.16) IV'ILp/2lvvi3p < c|v^|^^+^)|^|J+"(i-'^)-''. 

To gain something from this inequality we require aa + b < 1 which is equivalent to 

n n\ f n 2n \ n n n , , n 

+ 1 



b L/,|l-6 



3-- + ->a--:— +l-- + - = l-- + (l + a) . 
p qJ \q Sap J 3p q p q 

This in turn yields a < 2q/n which is certainly true by (|6.14p and ()6.15p . With 
K := aa + b ^ (0,1) we obtain the estimate 

|V$'(^(i))U,(o)" < C|^(0|^3(o)|^(t)|Si(n). 
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valid for a.e. t S [0,T] C [0,Tmax) and some m > 0. Similarly one obtains 

for a.e. t e [0,r] C [0,r,„ax)- Finally this yields 

(6.17) \<P'm))\HHn) < Cmtrm^in)mtKi(ny 

for a.e. t £ [0,T] C [0,Tmax)- Integration of the p-th power of (|6.17p and Holder's 

inequality imply ()6.12p . 

In conclusion we have the following result. 

Theorem 6.1. Let p>n + 2, n<3, q = max{2, j^g}, with 9 from (|6.5p . Sup- 
pose that a,c G C^(il)" and b G C^{fl) satisfy condition (H) as well as diva(x) — 
divc(x) = 0, X Cz fl and {a{x)\h'{x)) — {c{x)\i'{x)) = 0, x G dfl. Assume furthermore 
that $ e C'^~(R) satisfies (|6.4[) . (j6.5p and (I6.13p . Then there exists a unique global 
solution (jl'jf) of (IS.ip on Jq = [0,To], with 

and 

provided that the data are subject to the following conditions. 
(i) feLpiJo;Lpin))nL,{Jo;L2in)), 
(ii) geLpiJo;H^pin)), 

(iii) hi e Lp{Jo; Wp-^^''{T)) n LgiJo; i2(r)), 

(iv) h2 e ifpi(Jo; Lp(r)) n Lp( Jo; i^p"'/^(r)), 

(v) v^o e i3pV'^''(^). 

(vi) 9,,Vo = /i2|t=o- 
T/ie solution depends continuously on the given data and if f = g = hi — h2 — 0, 
the map ■00 '~^ V'(Oi ^ ^ ^+j defines a global semiflow on the natural phase manifold 
defined by (v) & (vi). 

Remark 6.2. The assertion of Theorem \6.1\ remains true if we assume that p > 
(n + 2)/3, which is sufficient for the well-posedness of the semilinear model by Remark 

7. Asymptotic Behavior 

In this last section we will give a qualitative analysis of global solutions of the Cahn- 
Hilliard-Gurtin system 

dtij} - div(a9tV) = div(W/x), t > 0, x^n, 

p.-c-Vn = I3dtip - A-ip + $'(-0), t>0, xen, 

(7.1) BVfi • z/ = 0, t>0, xeT, 

d^ip = o, t>o, xe r, 
0(0) = Vo, t = o, X en. 

To be more precise we will show that each trajectory converges to a stationary point, 
i.e. to a solution of the corresponding stationary system. The so called Lojasiewicz- 
Simon inequality will play an important role in the proof of this assertion. Assume 
that a,c e C^(17)" and b e C^(fi) with diva(a;) = divc(a;) = 0, x e il and 
{a{x)\i'{x)) — {c{x)\i^{x)) — 0, X £ dVt. Suppose that the data {(3,a,c,B) satisfy 
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condition (H) for all x G ft- Moreover we assume that $ G C'^(M) and that it satisfies 
the estimate 

(7.2) l*"'(s)| < C(l + |sP), for aU s e M, 

and some constant C > 0. Here 7 > 1 is arbitrary if n e {1, 2} and 7 < 3 if n = 3. 
At this point wc want to remark that (|7.2p already implies (|6.13p . 
Let tpo G ^IC^) s^'^h that di,^po = and let {ip,fJ.) be the unique global solution of 
(17.11) . We recall from Section IH] the energy functional 

E{u)^l\Vu\l+ f ^u)dx, 
^ Jn 

defined on the energy space 

V -.^{ue H^{n) : I udx^O}. 
Jn 

Note that due to (17. ip ^ and the boundary condition (I7.IP 3 we obtain J^tp dx = 

Jq ipo dx, since (a{x)\v{x)) = on F. If we perform a shift of ^ by means of ■0 = ip — c, 

where c :— J^^q dx, it follows that V' is again a solution of (17. ip . provided that the 

physical potential $ is replaced by $(s) = $(s + c). Additionally it holds that 

J^ip dx = 0. It follows from (|6.3p that E(tlj{-)) satisfies the equation 

j^Eiipit)) + m^it)\l + {a + c\^Mt)v^l{t)h + iBv^lit)\v^l{tyk = o, 

for all t £ IR+ . Making again use of Hypothesis (H) we obtain the inequality 

(7.3) j^Eim)+s{\dMt)\l + \V^,it)\l)<0, 

which holds for all t G R+. Integrating with respect to t and making use of (|6.4I) as 
well as of the Poincare inequality we obtain the a priori estimates 

V-e Loo(K+;i?l(f^)) and dt^p,\\7n\ e L2{R+ x ft). 

Proposition 7.1. The orbit {V'(0}teR+ ^s relatively compact in V. 

Proof. We rewrite equation (17. IL as follows 

/3dti'-AiP + ^lj = fi-Jl- (c(a;)|V/i) +7I+V-$'(0), 

where 'P = Tfi\ Jq *&'(V') dx. By the energy estimates above and the Poincare- Wirtinger 
inequality it holds that 

/ := ^ - 71+ (c|V/i) e L2(K+; L2{n)). 

Furthermore we have 

g := /!+ - $'(V^) e Loo{A+;Lg{n)), 

where q = 6/(7 + 2) is determined by the growth condition (|7.2p on $. The operator 
A := A — I with domain 

D{A) = {ue H^{n) ■.d^u = on F} 

generates an exponentially stable, analytic Co-semigroup {T{t)}teR^ in Lp{fl). There- 
fore 

T(.) */ e H^{R+;L2{n))nL2{R+;H^{n)) ^ Co{R+;H^m). 
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For the function g we apply elementary semigroup theory to obtain 

for each s E (0,2). The space H^{n,) embeds compactly into i?2(^)j if * i^ chosen 
close enough to 2. This completes the proof of relative compactness, since T{-)ipo e 

a 

The following proposition provides some properties of the w-limit set 

uj{tp) = {ip eV : 3 (tn) /• oo, s.t. i\j(tn) -^ ^ in V}. 

Proposition 7.2. Suppose that {tp, fi) is a global solution of (j7.ip and let $ satisfy 
Hypotheses (|6.4p and (|7.2p . Then the following statements hold. 

(i) The mapping t i— ^ E{il){t)) is nonincreasing and the limit limt_j.oo Eiipit)) =: 

Eca G K exists. 
(ii) The uj-limit setu){'ip) d V is nonempty, connected, compact and E is constant 

on uj{tp). 
(iii) Every ij^ao G ^(jp) is a strong solution (in the sense of L2) of the stationary 
problem 

(7.4) 

^ ' 9^Voo = 0, xeT, 

where fioo = mj In '^'(V'oo) dx = const. 
(iv) Each tpoo G '^(V') is o, critical point of E, i.e. i?'(V'oo) = m Y* , where V* 
is the topological dual space ofV. 

Proof. Inequality (I7.3p implies that E{'ip{-)) is nonincreasing with respect to t. Fur- 
thermore by (|6.4p it follows that E{u) is bounded from below for all u G V. This 
proves (i) . Assertion (ii) follows easily from well-known facts in the theory of dynam- 
ical systems. 

Let tAoo G "^(V")- Then there exists a sequence (i„) /^ -l-oo such that il^{tn) — >■ ^00 in 
V as n -^ 00. Since dtip G i2(R-(- x ^) it follows that ipitn -l-s) — >■ "000 in i2(^) for all 
s e [0, 1] and by relative compactness also in V. Integrating (|7.3p from i„ to t„ + 1 
we obtain 

E{i:{tn + I)) - E{i:{tn)) + e f f {\Vfi{tn + s,x)\^ + \dt^p{tn + s,x)\^) dxds<0. 

Jo Jn 

Letting t„ — > +cx} yields 

|VAi(t„ + •, ■)\,dtij{tn + •, •) ^ in L2([0, 1] x Ci). 

This in turn yields a subsequence (inj.) such that |V/u(inj. -|- s)|,9t'0(infc -I- s) — >■ in 
L2(^) for a.e. s G [0,1]. We fix such an s, say s* G [0,1]. The Poincare-Wirtinger 
inequality implies that 

Htn,+S*)-tl{tn,+S*)\2 
< Cp (\\/fi{tn, + S*) - VM(i„, + S*)\2 + J \<^'{'4'{tn, + «*)) - $'(VXi«, + S*))\ dx 

since J^fj. dx — J^^ <I''(V') dx. Letting k,l ^ 00 and making use of (|7.2I) it follows that 
A*(infc + s*) is a Cauchy sequence in L2(^), hence it admits a limit, which we denote 
by /ioo- Since the gradient is a closed operator in L2(^; R") it holds that ^00 G i?2 (^) 
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and V//00 — 0. Thus /ioo ~ const, and we have the identity /ioo = 701/0 '^'(V'oo) dx. 
Finally we multiply (jT.ip ^ by a function (p £ V in L2{^) to the result 

(7.5) {fl{tn, + S*),ip)2 + (C • V^(^„, + S*), ip)2 

= P{dMtn, + s*), ^)2 - (AV-Ci™, + s*), ^)2 + ('^'CV-Ci™, + s*)), ^h. 
Taking the limit t„j, — )■ 00 we obtain 

where a : F x F ^- K is defined by a{u,v) — (Vm, Vw)2 and (•,-)2 denotes the 
scalar product in L2(^)- Since $'(V'oo) € Lq{i^) with g = 6/(7 + 2) it follows that 
ipoc € D{Aq) = {u e Hy{V,) : d^u — 0}, where Ag is the part of the operator A in 
Lg{Q) which is induced by the form a{u,v). Observe that q > 6/5 by assumption, 
whence we may apply a bootstrap argument to conclude V'oo G ^2 (^) ^^id d„^poc — 
on r (recall that q > I may be arbitrarily large in case n E {1, 2}). Going back to 
(17.51) we obtain for (t„^ ) /^ 00 the identity 

(V-0oo,V^)2 + {<^' {l/j 00), (fi) 2 = (^^00,95)2, 

for all functions ip £ V. This yields (iii) after integration by parts. To prove (iv) 
observe that by I17[ Proposition 5.2] the functional E is twice continuously Frechet 
differentiable and its first derivative is given by 

{E'{u),h)v''.y ^ / WuWhdx+ / <i>'{u)hdx, u,heV. 
Jn Jn 

Integration by parts finally yields assertion (iv). 

D 

At this point we could simply refer to the paper of Miranville & Rougirel ^13 to 
prove the main Theorem 17.41 below. However, for the sake of completeness we provide 
a proof of this result. 

The next proposition is the key for the proof of the convergence of the orbit {Tlj{t)}t>o 
towards a stationary state as i — >■ 00. 

Proposition 7.3 (Lojasiewicz-Simon inequality). Let Lp e a;(V') and assume in addi- 
tion to (j6.4p and (17.21) that $ is real analytic. Then there exist constants s G (0, 2], 
C, (5 > such that 

\E{u)-E{p)\'-' <C\E\u)\v', 

whenever \u — if\v < 6. 

Proof. This is Proposition 6.6 in 0]. D 

Now we are in a position to state the main result of this section. 



Theorem 7.4. Let $ satisfy the conditions (|6.4|) and (17. 2p . Assume in addition that 
<& is real analytic. Then the limit 

lim 7/1 (i) =: V^oo 

i— fC30 

exists in V and ij^oo is a strong solution of the stationary problem (|7.4p . 



ON A CAHN-HILLIARD-GURTIN SYSTEM 33 

Proof. Since each element tp G a;('0) is a critical point of E, Proposition 17.31 implies 
that the Lojasiewicz- Simon inequality is valid in some neighborhood of (^ G ^{ip)- By 
Proposition [7?2](ii) the w-limit set is compact, hence there exists N gN such that 

N 

i=i 
where Bg {(fj) C V are open balls with center ipi G oj{4') and radius Si. Additionally 
in each ball the Lojasiewicz-Simon inequality is valid. It follows from Proposition [7?2] 
(i) and (ii) that the energy functional E is constant on Ld{tp), i.e. E{(p) = E^o, for all 
ip G u}{'ip). Thus there exists an open set U D U!{ip) and uniform constants s G (0, i] 
C,S >0 with 

\E{u)-Eoo\'-'<C\E'{u)\v', 

for all u G U. A well-known result in the theory of dynamical systems sates that the 
w-limit set is an attractor for the orbit {'0(i)}tgR^. To be precise this means 

lim dist(i/i(i),a;(V')) =0 in V". 

This implies that there exists some time t* > such that ip{t) G U for all t > t* and 
thus the Lojasiewicz-Simon inequality holds for the solution ip{t), i.e. 

(7.6) \EWt))-E^\^-'<C\E'(Mmv-*, t>t*. 

Define a function H : R+ ^ R+ hy H{t) = {E{tp{t) - E^f . Then with ([73]) and 
(frBj) it holds that 

^j^H{t) = {Emt)) - E^y-' (-jt^wt)) 



(7.7) > e 






The first Frechet derivative of i? in T^ reads 

{E\u),h)v*,v = / WuWhdx+ / <^'{u)hdx, 
Jn Jn 

for all (u, h) G V X V. Setting u = 4>{t) and making use of (|7.iP o we obtain with the 
help of Holder's inequality, Poincare's inequality and integration by parts 

{E'{ij{t)), h)v'.v = / (Kt) - Pit))h dx- c- \7fi{t)h dx - (3 dtipit)h dx 
(7.8) Jn Jn Jn 

<C{\V^l{t)\2 + \^Mt)\2)\h\2, 

since divc(x) =0, x G fl and {c{x)\iy{x)) = 0, a; G dfl. Taking the supremum in 
(|7.8p over all functions h G V with norm less than 1 it follows that 

\E'im)\v'<ci\v^iit)\2 + \^tiJit)\2). 

We insert this estimate into (|7.7p to obtain 

-^^Hit)>c,i\v^l{t)\2 + \^Mt)\2). 

Integrating this inequality from t* to oo it follows that \dti{j{-)\2, |V/^(-)|2 G Li(R-|_), 
since H{t) > 0. This implies that the limit limt_j.oo ipit) =: ipoo exists firstly in L2(^) 
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but by relative compactness also in V. Finally, by Proposition 17.21 (iii) the limit V'oo 
is a solution of the stationary problem (|7.4|) . The proof is complete. 

D 

8. Appendix 

Proof of Proposition \3.3[ 

We substitute p.9|) o into (|3.9|) ^ to obtain the elliptic problem 

(8.1) 

H+A{x,d)^i = div{adiv{DVu))-div{DVu)+f, a; G M"; B{x,d)^i = hi, x e dWl, 

with 

f = f3f + a-Vg-g&LpiJsxRl). 
Here the differential operators A{x, d) and B{x^ d) are defined by 

A{x, d)n := -(a + c) • V^ + div(a(c • V^t)) - div(/?BV/i), x e R!;:, 

and 

B{x,d)ii:=BV^i-v, a; G ^M" , 

valid for all ji G i7p(R"). It will be convenient to rewrite the operator A{x,d) as 
follows. A{x, d) — Ao[x, d) + Ai{x, 9), where 

Ao{x, d)n := - div [ l3BVfj. - -{a® c + c<g) a)V/^ j , x e R" , 

and 

Ai{x,d)n := a- {\7c\/^i) - - Div(a (g) c + c(g) a) • V^ - (a + c) • V//, x e R". 

Actually this splitting shows that problem (j8.ip is indeed elliptic by Assumption (H) 
and Proposition 12. 1[ provided w > is sufficiently small. Will will now proceed in 
several steps. 



Step 1. In this first step we want to reduce (|8.ip to the case of homogeneous boundary 
conditions B{x^ D)fi = 0. Consider the elliptic problem with constant coefficients 

(8.2) Xfi~div{BoVfi)^ f, xeRl, [BqV ^ji\u) = g , x e ^R" , 

where Bq := /3Bo — | (ao (8) cq + cq (8) ao) and A G R is a parameter. Note that (18.21) is 
an elliptic problem with a conormal boundary condition and constant coefficients. 
Thanks to Proposition 12 . 1 1 the matrix Bq is positive definite. By well known results it 
follows that for each / e Lp{Wl) and g G Wp~'^^''{dRl) problem (g^ has a unique 
solution fi G i/p(R" ), provided A > 0. We remind that the variable coefficients 

a{x) = oq + ai{x), c{x) = co + ci{x), B{x) = Bq + Bi{x), 

have a small deviation from the constant ones ao,ca, Bq, i.e. 

|ai|oo + |ci|oo + l^lloo < w, 

with oj > being sufficiently small. Furthermore we have a, c G VF^(R";R"), B G 
iyj^(R!f:;R"^"). Therefore we may apply perturbation theory to conclude that there 

exists Ao > such that for each / G Lp(R+) and g G Wp'^^^idWl) the elliptic 
problem 

(8.3) A^i - div(i3VAi) = /, a; G R", BV^i-u = g, x C, dW^, 
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has a unique solution /i e iJp(R"), provided A > Ao and cj > is sufficiently small. 
Here B := l3B - ^{a (g) c + c (g> a). Note that Ao{x, d)^i = - div(SV/i) and 

since (a\v) = {c\v) = by assumption. Moreover, the operator Ai{x, d) defined above 
contains only terms of lower order with Loo-coefficients. Thus, applying perturbation 
theory one more time, there exists Ai > such that for each / G Lp(R") and g £ 
Wl~^'^{dWl) the problem 

(8.4) \^i + A{x,d)n^ f, .tgM!^; B{x,d)fi^g, x E dR^, 

has a unique solution /i e iJp(R"), provided A > Ai and a; > is sufficiently small. 

Step 2. The results of this first step enable us to reduce (|8.ip to the case of homo- 
geneous boundary conditions. In this step we show that the Lp-realization A of the 
boundary value problem (A, B) with domain 

D{A) ^{ue i?p(M+) : B{x, d)u = 0}, 

is dissipative. First, let p > 2. Integration by parts yields 

Re / Aw wlwP"^ dx 



+ 
= - / \wY'~^ Re {j:{B\7w ■ \7w)\w\^ + (f " ^^^^ ' ^w)^"^) ^x 

for all w e D{A), since diva(x) = divc(x) = in 17 and {a{x)\v{x)) = {c{x)\v{x)) — 
on dVt. Setting Vw — u + iv and w = hi+ 162 with u,v E M", bj E M, we obtain the 
estimate 



Re (^{B\/w \/w)\w\^ + (^- l){B\/w\7w)w' 



> 



iiBu\u) + iBv\v))ibl+bl) 



>e(3{\u\^ + \v\^){bj+bl) = e(3\Ww\'\w\\ 

Here we made use of Proposition 12.11 This shows that A is dissipative in Lp{n) for 
p > 2. If p S (1,2) we replace \w\ by w^ :— -\/|it;P + e for e > in the calculations 
involving A and then pass to the limit as e \ 0. 

The dissipativity of A allows us to set A = 1 in (|8.4[) . By the same arguments one can 
show that the Lp-realization Aq of the elliptic boundary value problem {Ao,B) with 
domain D{Ao) — D{A) is dissipative, too. Indeed, Aqij, = div(i3V/i) and B defined 
above is a positive definite and symmetric matrix by Proposition 12.11 Therefore we 
may also set A = 1 in (|8.3p . 



Step 3. By the results of Steps 1 & 2 we may decompose the unique solution /^ E 
H'^-,{Wi^) of (|8.1[) into ^ = ^i + ^2, where ^1,^2 & ^p(R+) are the unique solutions 
of the elliptic problems 

(8.5) /ii+A(a;,a)/ii =div((adiv(DVu))), x eRI; B{x,^)^Il=0, x E dWl, 

and 

(8.6) 

l22 + A{x,d)iJ2 = f-diY{DVu)-Aiix,D)fXi, a; e ]R'|; B{x,d)ii2=g, x E dR\. 
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For fj,2 we have the estimate 

lM2U|(R:f) < C (^|w|i/2(R^) + |/|lp(R!; 



l5l, 



-i/p/ 



If^ilm 



with some constant C > 0, since Ai{x,d) consists solely of lower order terms with 
ioo-coefficients. To obtain the desired estimate for /i, we therefore have to prove the 
estimate 

for the solution /xi of (|8.5p . where C > 0. For this purpose note that the Lp-realization 
Aq of {Ao,B) generates a Co-semigroup in Eq := Lp(K") and A + Aq is a linear 
isomorphism from Ei := D{Ao) to Eq for each A S p{—Ao), the resolvent set of ^o- 
Let £^1/2 := [Ea,Ei]i/2, E_i/2 ■= (eI/^Y where S" = E' and denote by A_i/2 the 
i?_i/2-i'ealization of Aq. Here the symbol [•, •]i/2 denotes the complex interpolation 
functor of exponent 1/2. It follows from ^ Theorem V.2.1.3 & Corollary V.2.1.4] that 
the operator A_i/2 is the generator of a Co-semigroup with ^(^-1/2) == p(^o) and 
A + A_ii2 : i?i/2 —>■ E_i/2 is a linear isomorphism for each A G p(— Ao)- It remains 
to determine the spaces i?i/2 and E_i/2. To compute -E1/2, we have to interpolate 
Sobolev spaces involving boundary conditions. This has been done e.g. in [18j and 
[1] . Following these results it holds that 



-^1/2 — [Eo,Ei]i/2 = fl'p(R"). 

Actually in T this result was proven for C^-coefhcients but the result remains true 
for Vl^i,-coeflicients. From this characterization we obtain 



E-i/2 = {H'p'iK))' > - + - = 1, 1<P< 



1 



Set F = adiv(DVw) G H^iRl) and / = divF e Lp{Wl). Then fii e i/^(]R!f.) is 
a solution of the abstract equation fj,i + Aq/ii — f. We claim that this / can be 
identified with a linear functional in i?_i/2 (for which we will write / again). Indeed 
the mapping 



ip I— > 



fifidx^: (/,(^)£_i/2,Hi,, 



defines a linear functional on Hh{Wl), since by Holder's inequality it holds that 



\fU 



sup 

1,(1")^-' 



{fi^)E_u2,H\ 



= sup 

l'^lijl,(B")<l 



fifi dx 



<\f\i 



«+)• 



Integrating by parts we obtain furthermore 



|/b_i/2 = sup 

P + 

= sup 

Ivlffl (M")<1 

p + 



div Ftp dx 



Ivl, 



sup 



><i 



fifi dx 



sup 

p + 



(F|V(^) dx 



< 



sup |F| 



P + 



•)IvIhi,(R!?) - l-F^kpCR^iR"), 
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where we also made use of {F\i') = (a|i^) div(DVu) = 0. Since A_i/2 is the E_i/2- 
reahzation of Aq (hence an extension of Aq) with 1 G /?(— yl_i/2) = p{—Ao) and since 
/ = div(aAu) e E_i/2, we obtain a constant C > such that the estimate 

lMl|ifi(R!^) = IMiUi/2 < C'l/U-i/2 ^ C'|-P'Ilp(R!^;R") < C'IwIhKR!^), 

for the solution /ii € ifp(K" ) of (|8.5p is valid. From the estimates for ni and /i2 and 
the embedding Hp{M.'^) ^->- i?p(M"), we obtain a constant C > such that 

In the case that the functions depend on the parameter t it follows that the estimate 

(8.7) 

lMWIffi(R!^) < C [\fit)\L^(Rr^) + \9{t)\HHR1) + l^l(*)lvi/pi-i/P(K„-i) + |w(t)k2(R^)j , 

holds for a.e. t E J — [0,T] where the constant C > is uniform in t, since the 
coefficients of the differential operators considered above are independent of t as well. 
Taking the p-th power and integrating (j8.7p with respect to t, we obtain 

(8.8) 

lMlLp(J;ffi(R!f)) < C (^|/|lp(J;(R!^)) + |5|Lp(J;H1(R!^)) + 1^1 Ilp(,7;W^i-i/p(k„-1)) + l"lLp(J;if2(R^)) 

Finally the estimate for d^u in Lp{J; Lp(K" )) follows from p.lL and (|8.8p . The proof 
is complete. 
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